Diffusivity, mobility and time-of-flight measurements: 
scientific fictions and 1/f-noise realization 
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Influence of equilibrium thermal 1/f-type mobility fluctuations on time-of-flight measurements is 
considered. We show that it can explain experimental time dependencies of transient photocurrents. 
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I. INTRODUCTION 

It is well known that definite aspects of chaotic mo- 
tion of quasi-free charge carriers in many materials can 
be described in terms of diffusivity, D , and mobility, jj, . 
According to definitions of these quantities in statistical 
physics, if AX(t) is displacement, or path, of a carrier 
during time interval (0, t) , and t much exceeds charac- 
teristic memory time, tq , of carrier's random walk, then 



1 



(AX 2 (t)) = 2D 



t 



(AX(t)) f = fxf , 



(1) 

(2) 



where the angle brackets (...)o and (...)/ denote aver- 
aging over thermodynamically equilibrium statistical en- 
semble and non-equilibrium one in presence of sufficiently 
weak force / per carrier (to be precise, we consider pro- 
jections of vectors just onto direction of this force). 

Thus, D and fj, are averages over infinite variety of 
carriers. However, in practice many authors interpret 
relations (HJ-© in the form like 



AX(t) = lift + A d X(t) 
{A d X 2 (t)) f = 2Dt 



(3) 



where AdX(t) is "diffusive component" of the displace- 
ment. Thus, the ensemble averages D and \i are as- 
signed to a separate carrier as characteristics of its indi- 
vidual random trajectory! 

Such the arbitrariness, or "scientific fiction" , implies 
illusion of complete knowledge of non-equilibrium noise 
which stays the same as in equilibrium. Consequently, in 
view of the Einstein relation, D = T\x , one may think 
that at (j)(t) 3> 1 , where 



fVDi 
T 



the carriers' displacement is almost non-random, and 
therefore the time-of-flight measurements of transient 



photocurrents (see e.g. works [3-@ and references 
therein) should give nearly ideal responses to pulse-like 
or step-like excitation (see Fig.1 and Fig. 2). 

In reality one observes [l|-|5[ something qualitatively 
different, similar to what is shown in in Fig.l and Fig. 2. 

Such the difference is conventionally explained as result 
of random variations of time of flight (TOF) because of 
trapping and de-trapping of excess carriers by localized 
states (the so-called "dispersive transport", see e.g. 0, 
:6| and references therein). This explanation indeed is 
good in many characteristic cases. But in many other 
cases it seems unsatisfactory. Perhaps, by this reason, 
for instance, in Q the idea of "frozen" spatial variations 
of mobility was suggested. 

Quite another explanation arises when one understand 
that the fantastic expression © should be replaced by a 
realistic one, 



AX(t) = jl(t)ft + A d X(t) , 

( W) )/ = /*) 



(4) 
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where ju(t) is random quantity even in absence of the 
trapping, merely because a separate carrier has no defi- 
nite mobility (as well as a separate person does not have 
"mean life" and even may not know about it)! 

This statement is irrefutable, but its comprehension 
is prevented by prejudices of the naive probability- 
theoretical way of thinking instead of the fundamental 
statistical- mechanical one For the first this state- 
ment was argued in in the course of unprejudiced 
statistical analysis of random walks of atomic-size par- 
ticles and later in [l4| (or see [lj|) and in [l7l - [20| was 
confirmed, by example of fluid molecules, as consequence 
of the first-principle equations of statistical mechanics. 

Importantly, according to these investigations, a mag- 
nitude of fluctuations (uncertainty [To| ) of individual mo- 
bility of charge carrier is generally on order of it itself 
(i.e. ju — fx ~ fx). Therefore probability distributions 
of AX(t) taken at different values of the dimensionless 
force (j> = 4>(t) (see above) differ one from another not 
only by their mean positions but also by their widths 
and shapes, as it is shown in Fig. 3, so that real drift of 
a packet of even non-interacting carriers is accompanied 
by its spreading, approximately proportionally to both t 
and / . Particular results of such drift are the dark-curve 
time dependencies in Figs. 1-2. 
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FIG. 1: Ideal (fat grey curve) and actual (dark curve) tran- 
sient photocurrent responses to switching on illumination. 
The quantity r Q is is mean time of flight. 



The purpose of the present paper is just to demon- 
strate that these effects of the native carrier's mobility 
fluctuations resemble effects of the trapping. 



II. WHY AND HOW MOBILITY FLUCTUATES 
BY ITSELF 

Introducing related fluctuating diffusivity D(t) = 
TJl(t) , one can ask equivalently: why carrier's diffusivity 
fluctuates by itself? The cause was explained in the men- 
tioned works. It is temporal non-locality of such physical 
characteristics as diffusivity: the latter by its sense de- 
scribes random motion of a carrier (or a particle in gen- 
eral) during time intervals significantly greater than the 
memory time r . Therefore at the end of any such in- 
terval the system already does not remember what took 
place at its beginning. Consequently, the system can not 
keep a definite value of the diffusivity. 

If it is seen, then it is clear that resulting fluctuations 
(uncertainty) of diffusivity (and hence mobility) , firstly, 
by their order of magnitude are comparable with its av- 
erage value. Secondly, they can be arbitrary long, that 
is have no characteristic time scale, and thus represent a 
"flicker noise" whose power spectrum Sd{ijj) unbound- 
edly grows at u) — > . Thirdly, their probability distri- 
bution is close to a stable one. 

Additional reasonings can be found in [IcHl^ , IT3 - 
[24[. To confirm them at fundamental level, we 
have to calculate statistical characteristics of the dis- 
placement AX(t) , such as its characteristic function 
(exp [ikAX(t)]} , merely following rigorous recipes and 
equations of the statistical mechanics. Or, at least, ana- 
lyze admissible forms of (exp [ikAX(t)]) at phenomeno- 



FIG. 2: Ideal (fat grey curve) and actual (dark curve) tran- 
sient photocurrent responses to short pulse of illumination. 
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FIG. 3: Comparison of evolutions of the AX(t) distribution, 
considered as a function of the dimensionless parameter (j> , 
in the conventional "fantastic" model Q (four identical sym- 
metric bell-shaped curves) and in realistic model Q accord- 
ing to results of [Tol. [Til. (dark asymmetrical curves). The 
most left bell-shaped curve corresponding to = is com- 
mon for both models. The two most right fat curves demon- 
strate qualitative difference of the models at <f> = 12 . 



logically under minimum of assumptions. 

For the first such investigation was made in 
(or see [l6j]) where thermodynamically equilibrium ran- 
dom walk AX(t) of an electron (or other particle) in 
spatially uniform medium was considered basing on the 
only assumption that it possesses diffusive scale invari- 
ance, AX(st) ~ y/sAX(t) (the tilde symbolizes statisti- 
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D/D = ji/fi 



FIG. 4: The three dark curves represents the diffusivity and 
mobility probability density distributions following from Eq[7] 
[II], GJ, El at different values of In (t/r ) and rl/Dro . The 
fat grey curve illustrates the diffusivity and mobility distri- 
bution (10) [17H20|. All three distributions have power- law 
long tail with the exponent between 2 and 3 (somehow cut off 
far on the right). The inset: the same in double- logarithmic 
scale. 



cal equivalence) but destroyed at some microscopic time 
scales To and tq (e.g. mean free path or mean hop length 
of the particle under consideration). It was shown that 
at t ^> T and \k\r <C 1 the characteristic function of 
AX(t) can be approximately expressed by 



(exp [ikAX(t)]) = exp 



-DkH 



In {no ft 



]k 2 ) 



In (r /i) 



(5) 



The corresponding statistics can be named "quasi- 
Gaussian" QlQji. 

Similar analysis of non-equilibrium situation, when the 
particle is affected by an external force / , needs in the 
help of the "generalized fluctuation-dissipation relations" 
[H HI] which show [HH3 H El that at sufficiently 
weak force, in the sense of fro/T -C 1 , the characteristic 
function can be obtained by simple replacement in ([3]), 
namely, 



-k 2 + ikf/T 



(6) 



(thus the Einstein relation extends from mean values to 
fluctuations of diffusivity and mobility). 

Therefore, at <f>(t) 3> 1 , neglecting in (TJ| the diffusive 
part AdX(t) of AX(t) , one can write approximately 



(exp [ikAX(t)]) = exp 



ik/ift 



\n(T /t-ikfr 2 Q /T) 
In (t /£) 



(7) 



The corresponding statistics can be named "quasi- 
Cauchy" since the corresponding probability distribu- 
tions (see dark curves in Fig. 4) are similar to the asym- 
metric stable Cauchy distribution [27], HU (see Fig. 5). 



FIG. 5: The maximally asymmetric stable Cauchy distribu- 
tion. It is principal ancestor of the distributions shown in 
Figs. 3 and 4. 



Just such statistics implies the displacement distribu- 
tions shown in Fi g. 3 b y dark curves. Its physical meaning 
was 

explained in [HQS 0. 
One can see that characteristic functions like ([5]) and 

([7]) mean presence of "flicker", or 1/f- noise type, fluc- 
tuations in diffusivity and mobility, with low-frequency 
power spectrum |9rll2l [la ] 



Sd 
D 2 



a(oj) 



2r 2 



3Dtq In (tqw) 



(8) 



This theoretical estimate very well agreed with experi- 
mental data in 1983 [TJ| and equally well agrees with new 
data about semiconductors [2!|, carbon nano-tubes [3fj| . 
graphene layers and graphene transistors [3ll [32[ and 
other electronic devices (at that revealing origin of typical 
values of the "Hooge constant", a(u) ~ 0.002 [l|]). It 
should be underlined again and again that this 1 /f- noise 
is by its nature an integral part of the atomic thermal 
motion and exists even in ideal crystal structure in ab- 
sence any traps (see explanations in [To|. [T2I. [IFjL [2ll. |22| ) ! 

The fundamental microscopic approach to analogous 
diffusivity and mobility fluctuations so far was tried only 
in case of molecular random walks in fluids Il4l420l |. on 
the base of the exact Bogolyubov equations [33| . 

In [3 (see also [HGJ]) it was found, for atoms of a 
gas, that 



S d {uj) Sp(u)) 



D 2 



UJ TqUJ 



(9) 



which satisfactorily agrees with experimental data on 
electrolytes [TH, [34[ (where the Hooge constant is anoma- 
lously large). 
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FIG. 7: Transient photocurrent response to step pulse illu- 
mination corresponding to the mobility distribution (10) at 
different /3 values. In the inset: the same in double logarith- 
mic scale. 



But the underlying probability distribution of diffu- 
sivity and mobility remained unknown. It was found in 
[T7L [20j by means of more rough but deeper penetrating 
approximation, and its density can be expressed as 



Up(x) 



x P+»T(fi + l) GXP I X 

b 



(10) 



(3D fin ' 



with /3 > . The spectrum © arose in the case when 
(3 = 1 [17]. The corresponding probability density U\(x) 
is shown in Fig. 4 by far grey curve. Evident property 
of all obtained distributions is that their most probable 
values are significantly smaller (three times at (3=1) 
than their mean values. 

Thus, both phenomenological and microscopical ap- 
proaches to the native thermal mobility fluctuations 
(when being correctly performed) gave qualitatively iden- 
tical results. 



III. MANIFESTATION OF THE NATIVE 1/F 
MOBILITY FLUCTUATIONS IN TOF 
EXPERIMENTS 

The TOF measurements of excess (injected) probe par- 
ticles, e.g. photo-excited charge carriers, usually presume 
that 4>(t) ^> 1 . According to the aforesaid, this means 
that S^f 2 > 2Tfi , that is at corresponding time intervals 
contribution to particle's path fluctuations from 1/f dif- 
fusivity/mobility fluctuations exceeds contribution from 
diffusion itself (i.e. from white thermal noise). In other 
words, TOF experiments are first of all observations of 
the native diffusivity/mobility 1/f-noise. 

Hence, it is reasonable first to consider such imaginary 
(but not fantastic) situations when the above mentioned 
simple theory is sufficiently adequate. This would deter- 
mine a background of the TOF spectra which exists even 
in absence of specific effects like the "dispersive trans- 
port" (transport via a large number of free electron traps) 
jot taken into account by the theory. 

At 4>(t) ^> 1 , if knowing the mobility probability dis- 
tribution, one can easy find that of TOF r = L/jif 
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(where L is distance to be overcame) and resulting time 
dependence of transient current (flow), and vice versa. 
The distributions corresponding to characteristic func- 
tion Q have no simple analytic expression but can be 
approximated by simple related distributions (10). Then 
calculation of of transient current become so trivial that 
its description can be omitted. 

The results for short-pulse and step-like excitations 
are demonstrated on Fig. 6 and Fig. 7 where the currents 
are normalized to unit start value and unit steady-state 
value, respectively (t/ there is mean TOF value). 

IV. DISCUSSION AND CONCLUSION 

Visual comparison of our Figs. 6 and 7 with experimen- 
tal data graphically resented in [H-Q (and similar works) 
shows rather high degree of correlation between one and 
another. This, gives evidence that many observations 
can be explained by equilibrium mobility/diffusivity 1/f- 
type fluctuation accompanying usual, "non-dispersive" , 



transport. On the other hand, we can say that many 
observations, when even their authors have doubts as to 
"dispersive transport", do proof that the latter is not 
unique cause of TOF and mobility distributions. 

For instance, may be.the experimental mobility dis- 
tribution found in work [2| (which stimulated the present 
work) should be assigned not to spatial inhomogeneities 
but to temporal mobility fluctuations of charge carriers 
(otherwise, why this distribution so strongly resembles 
the picture in Fig. 4 ?). Analogously, the large width and 
asymmetry of mobility distributions of ions in air mea- 
sured in [35j j can be manifestation of flicker fluctuations 
in ions' mobilities. 

At the same time, undoubtedly, our theoretical gallery 
does not include specific pictures of pronounced "disper- 
sive transport". This means merely that 1/f- noise under 
dispersive transport has more rich scaling properties and 
hence more complicated mobility statistics than what 
was considered; Developing of theory of such 1/f-noise 
from rigorous statistical mechanics, instead of probabilis- 
tic "scientific fiction" , is interesting task for the future. 
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